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Abstract 

We study the Toda equations in the continuous level, discrete level and ultradiscrete 
level in terms of elliptic and hyperelliptic a and if) functions of genera one and two. 
The ultradiscrete Toda equation appears as a discrete- valuation of recursion relations 
of ip functions. 

1 Introduction 

Recently Kimijima and Tokihiro found solutions of the discrete and ultradiscrete Toda 
equations in terms of elliptic and hyperelliptic functions in this article we present 
another type of solution of another type of the discretization of the Toda equation [7j and 
its ultradiscretization ^5] associated with algebraic curves of genera one and two. 

Elliptic and hyperelliptic a functions are related to nonlinear differential equations from 
the beginning We study the Toda equations at the continuous level, discrete 

level and ultradiscrete level in terms of a functions and ip functions of genera one and two. 
We show that these equations have solutions expressed in terms of a and ifi functions. 
Here the ifi functions are defined by rational functions of the a functions. 

In |LSj it was shown that the ip functions can be related to discrete nonlinear equations, 
such as the discrete Painleve equations. This article can be considered as one of a series 
in which relations between ip functions and discrete nonlinear difference equations are 
unfolded. 

Further, as was mentioned in ^3], the ultradiscrete sometimes can be regarded as a 
valuation of a related field. This article shows that, in the case of the Toda equation, it 
can be also realized as a discrete valuation of a function field over a Jacobi variety. 

In Section 2 we concentrate on the genus one case and give concrete solutions. We 
investigate the genus two version in Section 3. It is shown that all solutions of the Toda 
equations in this study are connected with the addition formulae of the a functions. 
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2 Genus One Case 

In this Section we deal with an elliptic curve given by 



Ci : \y 2 = y 2 = x 3 + X 2 x 2 + Xix + A 



= (x-h)(x-b 2 )(x-b 3 ), (2.1) 
where the b s are complex numbers. 

2.1 Continuous Toda Equation 

Firstly we give a p function solution of the continuous Toda equation ^H] . We treat the 
Weierstrass elliptic a function associated with the curve C\, which is connected with the 
Weierstrass p function by 

p(u) = -^log ff (u). (2.2) 
A local parameter u in C\ is given by 



( x '») dx 

and p(u) is equal to x(it). Here oo is the infinity point of C\. 
The addition formula of the a functions is given by 



(2.3) 



a(v + u)a(u — v) 
[a(v)a{u)f 



- p(u)] = FTT^bb^ ' ( 2 - 4 ) 



By differentiating the logarithm of (|2.4j) by n twice, we have 
d 2 

- j"2 log[p(u) - p(w)] = p(« + v) - 2p(n) + p(n - «). (2.5) 

For a constant number Mo, by letting u = uuq + 1, v = uq and b := p(uo), we have 
d 2 

~dt 2 log [ p ( nu ° + ~ 6 1 = [^(( n + 1 )'" + *) ~ h l 

-2[p(nn + t) - b] + [p((n - l)n + t) - b]. (2.6) 

Further by letting q n := log[p((n + 1)uq + t) — b], we have 

d 2 

= e9n+1 ~ 2e9n + e9n ~ 1 - (2 - 7) 

This is identified with the continuous Toda lattice equation. In fact, by letting q n = 
Qn — Qn-i, Qn obeys the nonlinear differential equation of a nonlinear lattice |16| . It is 
clear that this elliptic solution comes from the addition formula (|2.4|) . Later we show that 
a genus two solution of the Toda equation can be expressed by a similar form in §3-1. 
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2.2 Discrete Toda Equation and ip Functions 



Though there are several models of the discrete Toda equations, we concentrate on a model 
given in 0. In this subsection we give elliptic solutions of the discrete Toda equation. 
The elliptic ifi function is given by 



a(nu) 
a(u) n2 



and, due to the addition formula (|2.4|) . it satisfies the recursion relation [T8| . 

m—n 



(2i 



(2.9) 



Further ip n can also be computed using the Brioches-Kiepert relation 0111) 



\n-l 



[l!2!---(n-l)!] 2 



p'{u) p"(u) 
p"(u) p"'(u) 



p 



(n-1). 



P (n \u) 



P [n \u) 



p 



(2n-3) 



(U) 



(2.10) 



Here derivatives of u are denoted by ' and <yn \ By noting 4- = 2y^, the ip function is a 
polynomial of x and y over the complex number C. In fact the ip function can be explicitly 
obtained as, 

^i(n) = 1, 
i) 2 {u) = -2y, 

^ 3 (ii) = 3x 4 + 4A 2 x 3 + 6Aix 2 + 12A x - Ai 2 + 4A 2 A , 
^ 4 (n) = -4y[x 6 + 2A 2 x 5 + 5Aix 4 + 20A 2; 3 

+(20A 2 A - 5Ai 2 )x 2 + (8A 2 2 A - 2A 2 A X 2 - 4AiA )x 

+4A 2 AiA -Ai 3 -8A 2 ]. 



(2.11) 



For t/j n , n > 4, we have the recursion relations 

^ 2n+ l = Vn+2V>n - V'n+lV'n-l, 

Thus we know that 

ip n (u) 6 C[x, Ao, Ai, A 2 ] for odd n, 
ip n (u) € C[x, Aq, Ai, A 2 ]y for even n. 



(2.12) 



(2.13) 



When n is odd, tp n is a polynomial of x whose order is (n 2 — l)/2 and, for an even n, the 
order of x for ip n /y is (n 2 — 4)/2. For specific curves, we give explicit forms of i/) n in the 
Appendix. 

We comment on the properties of tp functions. We note that cx(u) is characterized by 
the property that it has no singularity with respect tow£C and its zeros are identified 
with a lattice points generated by the periodicity 2u> of p(u), p(u + 2w) = p(u). In 
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other words the zero of a is congruent to the origin of the local parameter u modulo the 
lattice. Accordingly, as tp function is a function over the curve C, we conclude that a point 
satisfying 

Mu) = (2.14) 

is a point for which nu is equal to the lattice point again. In other words we have n-cyclic 
points as zeros of ip n . Conversely it can be shown that the polynomial of x and y whose 
zeros multiplied by n are lattice points must be ip n modulo constant factors. 
Hence, if n is a factor of m, i.e., n\m, it is clear that ip n is divided by ip m , 

Mm- (2.15) 
For mutually coprime numbers p, q and an integer no, we introduce 

$i '■= ipno+pi+qj- (2-16) 
By letting n = uq + pi + qj we have 

1p n+p 1p n - p = 1pllp p+l 1p p -i - ijjptpn+ilpn-!, 

Ipn+qlpn-q = ^n^q+l 4>g-l ~ 4>q • (2.17) 

The components ipn+fpn—i m both formulae give a relation, viz 

(V'n+pV'n-p - ifyp+l^p-l)^ = (i>n+q1pn-q ~ ifi^q+X^q-l)^. (2.18) 

Noting n = hq + pi + qj (|2,18|) can be regarded as an evolution equation for i and j 
when we consider ip q , ip p and ip q ±i as initial values. Let us assume that ip q ip p and 
ipp± q are not equal to zero by choosing the parameter u. By letting 5 := ipq/ipp and 
c(l - 5) := -ipp +q ipp-q/(ip q ) 2 we have 

" < l " #)<l>Hi S ~ * 2 ^M-i = 0. (2-19) 
For later convenience we do not fix c and define 

Vj:=(^f]-c (2.20) 



Then we obtain 



log | <° + V ' >' , (C+ .^' ) ' ■ |. (2.21) 



When c = 1, this equation is one of discrete versions of the Toda equation, which appeared 
in 0. 

The condition c = 1 means that 



Ipp+qlpp-q + 4>p1pp ~ Ipqlpq = 0, 



(2.22) 
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which is an equation of (p 2 + q 2 — 2)-order with respect to x. In other words, for a point 
u satisfying (|2.22|) . we have a solution of the discrete Toda equation in [7] in terms of tp 
functions. 

Further we introduce U, 3 := (V, 1 + c) which satisfies 

( W? \ = ( (cp^-fHfu^l 

Wi^Wt 1 ) ) \«l - 5 2 ) + PUilJWL - 5 2 ) + 

We investigate this equation with general c and go on to call it the discrete Toda equation 
in this article. We note that this solution is due to the recursion relation ()2.9|) which comes 
from the addition formula l|2.4[l . 

2.3 Periodic Solutions of Discrete Toda Equation 

For general c in (|2.2(Jf) we consider a periodic solution of (|2.23() . It is obvious that, when 
ip n = 0, ipnr = 0. Thus there may exist a point, u\, such that 

i>i{u 1 ) = il> n + i {u 1 ). (2.24) 

In fact we have solutions of (|2.23|) for a curve y 2 = x 2 (x + 1/4) and its point x = —1 (See 
the Appendix). The ip function has values as in Table 1. 

Table 1: i\) n at x = — 1 



01 2 3 4 567 8 9 10 11 12 





For (p,q) = (3,2) and no = 0, i.e., 5 2 = —3/4 and c(l — 5 2 ) = 1/4, we have a periodic 
solution of (E2SI). 



Table 2: Tj{ {p = 3, q = 2 case) 



j\i 





1 


2 


3 





oo 





oo 





1 


1/3 


3 


1/3 


3 


2 


1/3 


3 


1/3 


3 


3 


oo 





oo 






For (p, q) = (2,3) and uq = 0, i.e., 5 2 = —4/3 and c(l — <5 2 ) = 1/3, another periodic 
solution of (|2.23l) is given in Table 3. 



Table 3 


: [// (p 


= 2, 


g = 3 


case) 


J'V 





1 


2 


3 





oo 








oo 


1 


1/4 


-2 


-2 


1/4 


2 


oo 








oo 


3 


1/4 


-2 


-2 


1/4 
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2.4 Ultradiscrete Toda Equations 

In this subsection we investigate the ultradiscrete version of the Toda equation using ip 
functions. 

For the elliptic curve C\ a local parameter t should be characterized by 



for a generic point xq in C±, 

for a finite branch point b% in C\, 

for the infinity point oo in C\, 



t = X — Xq, 

t = y/x~^b h (2.25) 

t = 1/y/x. 



Let a localization of the commutative ring R = C[x,y]/(y 2 — f(x)) at u be denoted 
by R uo . Let K uo be a field of Laurent transformations at uq of rational functions. The 
valuation of the field K Uo is given that for / G K Uo , let val(/) = oo if / = 0, if / is given 
by 

f(u) = a{u - u ) m + 0{(u - n ) m+1 ) (2.26) 

for a / 0, let val(/) = m [£]. Denoting set of integers by Z, the discrete valuation is 
known as a map 

val : K uo -> Z + oo, (2.27) 
which satisfies 

val(/#) = val(/) +val(g), 
val(/ + 5 ) > min(val(/),val( 5 )). (2.28) 

For example the inequality in (|2.28j) appears due to a case, k = m and a = — b for 
/ = a(u — uq)" 1 + • • • and g = b{u — uo) h + • • • with (a, b ^ 0). Inversely, as long as we 



avoid such a case, we can regard the second relation in (|2.28|) as an equality. 
R Uo can be expressed as 

R uo = {f£K U0 \ val(/) > }. (2.29) 

Ruo := { / ^ K-ua I va K/) = } is a multiplication group in R UQ . An element in R* is 
called unit. The subset m := { / G K UQ \ val(/) > } of R UQ is a unique maximal ideal 
in R Uo and thus we have a filter structure, 

m k D m fe+i_ ^ 30 ^ 

Here the multiplication among ideals is given as a set of sum of multiplications of elements 
in the ideals. Due to the filter structure there naturally appears a nonarchimedean distance 
given by 

|/ - <7| val := exp(-val(/ - <?)). (2.31) 

Thus an element / in m is a smaller element than unity, i.e., |/| va j < 1- When 5 behaves 
like a small parameter 0, we regard it as an element of m, i.e., 

5=^%)Gm. (2.32) 
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We now consider the point satisfying 



c(l - 5 2 ) 



(2.33) 



Define 



fi j := -val(C^), d := -val(<5 2 ). 



(2.34) 



When we expand them as U-5 2 = a(u — uo) m H and c(l — 5 2 ) = 6(n — uo) fc + • • • with 

a,b 7^ 0, we assume that for any i and j, A; is not equal to m or a is not equal to —b if 
k = m. Then (|2.23j) becomes 

- 2/ 4 J + f^ 1 = max(0, /^^ + d) - 2max(0, / 4 j + d) + max(0, f.^ + d). (2.35) 

This is identified with the ultradiscrete Toda equation in |15| . 

Let us consider solutions of the ultradiscrete Toda equation (|2.35|) . By letting g n := 

val(VvO, 

For the curve y 2 = x 3 + 1/4, and uq at x(uq) = (— 1/4) 1 / 3 , we have gi as in Table 4: 
Table 4: g n at y = 



n 


12 


3 4 


5 


6 


7 8 


9 10 11 12 ••• 




oo 1 


1 





1 


1 


1 1 


Then we have a 


solution of (|2.35|) for p 


= 3, 


9 = 


= 2 and 


n = 0: d = -2, val(c(l - 5 2 




Table 5: // (p = 


= 3, < 


7 = 


2 case) 








1 


2 


3 


4 


5 ••• 







CXD 


-2 


2 


-2 


2 ••• 




1 


2 


-2 


2 


-2 


2 ••• 




2 


2 


-2 


2 


-2 


2 ••• 




3 


2 


-2 


2 


-2 


2 ••• 



Next we deal with a curve y 2 
$iven in Table 6. 

Table 6: g n at x = 



x and a point uq(x = 0). The values of gi are 



n 



0123456 7 8 9 10 11 12 



9n 



oo 1 4 5 8 13 16 21 28 33 40 49 
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When (p, q, no) 



(5, 2, 0), we have val(c(l - 5 2 )) = 0, d = 14 and Table 7. 
Table 7: f- (p = 5, q = 2 case) 



j'V 


1 


2 


3 


4 





oo 


18 


14 


18 


i 


18 


14 


18 


18 


2 


14 


18 


18 


14 



In this case, |<5| va j > 1. 

3 Genus Two Case 

In this section we investigate genus two solutions of the Toda equations using the hyper- 
elliptic cr functions and ip functions. 

The hyperelliptic cr function was defined by Klein after the prototype had been discov- 
ered by Weierstrass ■ Let us fix a hyperelliptic curve with genus two, 



C 2 : y 2 = x 5 + X 4 x 4 + X 3 x 3 + X 2 x 2 + X lX + X , 



(3.1) 



where Aj, i = 0, 1, ••• ,4 are complex numbers. For a point in the symmetric product 
space of the curve C 2 , ({x\, yi), (x 2 , 2/2)) £ Sym 2 C2, its corresponding point u = (ui,u 2 ) 
in the Jacobi variety J 2 is given by 



'(xim) dx f^ y2) dx 
Uf.= I VI — , u 2 :-- 

y Joo y j 00 



{xi ' y ^ xdx f ix2 ' y2) xdx 

y Joo y 



(3.2) 



Here 00 means the infinity point of the curve C 2 . At the point, p functions of genus two 
are defined as 



P11 



f(xi,x 2 ) - 2yvy 2 



{xi - x 2 ) 

<2 



Pl 2 = XlX 2 , p 22 = xi+x 2 , 



(3.3) 



where f(x, z) := X^j=o xJ zJ (^2j+i{ x + z) + 2X 2 j). It is known that there is a global function 
over C 2 such that 



pij 



& 2 



log cr, 



(3.4) 



duiduj 

which is the a function of genus two. 
3.1 Continuous Toda Equation 

Though there were found solutions of the continuous Toda equation in terms of the 9 
function related to a hyperelliptic curve of genus g in |JJJ, in this article we give another 
type of expression of solutions in terms of a functions related to a curve with genus two. 
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The additive formula of a function of genus two is given by [2] , 
a(v + u)a{v - u) 

By letting 

Q(u,v) := -(pn(u) - pu(v) + pi 2 (u)p22(w) - Pl2(v)p22(u)), (3.6) 
we have 

() 2 



duiduj 



log(Q(u,v)) = pij(u + v) -2pij(u) + pij(u-v). (3.7) 



Let us fix u = uuq + t, v = uo, constant numbers bij := Pij(uo) and 
d 2 d 2 d 2 

A := ami + b22 m^r 2 + bl2 d^df 2 (3 - 8) 

Then we have a relation, 

Q((n + l)uo + t,uo) — 2Q(nuo + t,u Q ) + Q((n — l)uo + t,uo) = -AlogQ(nu + t,u ). (3.9) 

By considering the relations (|3.3|) we let uq correspond to a point ((a?i, yi), (x2, ^2)) £ 
Sym 2 C2 and then have 

^22=^1+^2, h 2 = X\X2- (3.10) 

If the points are mutually conjugate or identical, i.e., x\ = X2, 

Hence for t := t\ + t2/xi, 

q n :=logQ(nc + t,c), (3.12) 
obeys the continuous Toda equation, 
d 2 

— —q n = e q " +1 - 2e qn + e 9 "" 1 . (3.13) 
dt z 

As we showed in the genus one case, this genus two solution also comes from the addition 
formula ljH.5JI . 

3.2 Discrete Toda Equation 

We give relations between the discrete Toda equation and tp functions of genus two as 
follows. Generalizations of the ip function in ()2.8|) to genus two curves are given by two 
different definitions; one is defined over the Jacobi variety J2 and another is defined over 
the curve C2. The former is studied by Kanayama 8 j and the latter is investigated by 
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Grant, Cantor, Onishi and this author (see the references in ^SJ). The definition by 
Kanayama is 8 

Mu) = 3=£. (3.14) 
a(u) n 

Further he showed that ij) n obeys the same recursion relation as (|2.9|) . viz 



m—n 



(3.15) 



basically using the the addition formula 1)3.5)1 . Hence tpk obeys a relation which has the 
same form as ()2.12p , Kanayama gave the explicit forms of ipi, ip2, ^3 an d ipi in terms of 
p functions (|3.3|) in [S] . We can compute an explicit form of any ip n as a rational function 
of the affine coordinates (£1,2/2) and (£2,2/2) of the curves Sym 2 C2 even though it is too 
large to give its explicit form here. 

Due to its form, it is obvious that 1)3.15)1 is also related to the discrete Toda equation. 
For mutually prime integers p, q and an integer no, we define quantities, 

0/ := ip no+pi+qj , (3.16) 

5 := ipq/t/jp and c(l — 5 2 ) = ipp+gipp-g/ipp 1 . Then 1)3.15)1 becomes 

Hence we have a solution of the discrete Toda equation ()2.21)1 in ^7 as shown in Section 
2.2. 

As a simple Abel variety has a division field as its endomorphism in the category of the 
Abel variety as it is known as Poincare's complete reducibility theorem Hence, even 
though the Jacobi variety J2 is two-dimensional, an isometry ip : J2 — > Ji is characterized 
by an integer. The zeros of ip n belonging to Z/nZ determine the isometry. Thus, as long 
as we deal with isometries of Jacobi variety, an extension of the ip n functions to functions 
with double-index must fail. It implies that ()3.16)) is a natural in the sense of a realization 
of the discrete equation in category of the Abel variety. 



3.3 Ultradiscrete Toda Equation 

We consider the Jacobi variety J2 as a commutative ring and its localization ring R UQ 

and/or a field K UQ related to R UQ . Similar to the case of genus one, we deal with a point 
of curve satisfying 

c(1 _ j2) = ^feH e ^ (318) 

By letting 

/ t J :=-vai r' — p j, d:=-val(5 2 ), (3.19) 
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and being supposed that for all of i and j, the valuations of the additions of / s expressed 
by the minimal functions like the equality case in the second relation of (|2.28|) . we find a 
solution of the ultradiscrete Toda equation in |15j . 

fj +1 - 2f t j + f^ 1 = max(0, + d) - 2max(0, f t j + d) + max(0, /.^ + d). (3.20) 

Even though the case of genus one gives us trivial solutions, genus two case is expected 
to provide us nontrivial solutions because it has larger degree of freedom than the elliptic 
curve case. 



4 Discussion 

In this article we have considered the relations between the Toda equations in the con- 
tinuous, discrete and ultradiscrete levels and a functions of genera one and two. We 
showed that these solutions, in principle, come from the addition formulae of the algebraic 
functions over algebraic curves of genera one and two. 

As we started from the curves, all of the solutions are expressed by points at curves 
(|2.1|) and IjH.lJ) . As we gave some explicit solutions related to elliptic curves, we can 
basically find explicit forms of the other solutions in terms of points of the curves even 
of genus two, though they might be slightly complicated. As a next step, we should give 
more explicit computations of the ip functions on the genus two case by finding a nicer 
strategy to handle the huge polynomials. However, it is remarkable that in our solutions, 
there do not appear excess parameters except the coefficients Aj, i = 0, 1, • • • ,4 in (J2.1|) 
and (j.S.lj) . In other words we have no ambiguity for the parameters even in genus two 
case and it means that we are free from the so-called Schottky problem. This contrasts 
with the solutions in terms of the 9 functions over the Jacobi variety El • Of course 
as it might be difficult to deal with the hyperelliptic integrals, thus we should select the 
methods according to the circumstances. 

Further it is interesting that the ultradiscrete equations can be defined on the Jacobi 
varieties associated with nondegenerate algebraic curves over a field with character zero 
using the concept of discrete valuation. (In ^3] we show that the ultradiscrete equations 
can be defined over fields with nonvanishing character.) It means that, if we find a recursion 
relation over an algebraic variety, we might have its ultradiscrete version by taking its 
discrete valuation. 

Finally we comment on the higher genus case. Unfortunately since the addition formula 
is not simple ^ 0], we could not deal with a functions associated with a curve with a 
higher genus as mentioned above. We hope that we can obtain such solutions in future. 
We note that the paper 0] may have some effects on the study. 

Acknowledgement 

I am grateful to Prof. Tokihiro for sending me his interesting paper with Kimijima |17| . 
I thank to Prof. Onishi for directing me to the theory of Poincare in |llj . 



12 



S. Matsutani 



Appendix 

Let us deal with y 2 = x 3 + 1/4, y 2 = x 3 — x and y 2 = x 2 {x + 1/4) and show explicit 
function forms of their ip functions. 



A.l y 2 = x 3 + 1/4 



^ = 1, (A.l) 

V>2 = -2y, (A.2) 

^ 3 = 3x(l + x)(l-x + x 2 ) =3x(l + x 3 ), (A.3) 

V>4 = M-^ + Wx 3 + 2x 6 ), (A.4) 

V> 5 = -l- 25x 3 - 15x 6 + 95x 9 + 5x 12 , (A.5) 
V>6 = V2^ 3 (-2 + x 3 )(l-3x + 3x 2 + x 3 ) 

x (1 + 3x + 6x 2 + llx 3 + 12x 4 - 3x 5 + x 6 ), (A.6) 
V> 7 = (l-x 3 + 7x 6 ) 

x (1 - 48x 3 - 741x 6 - 1924x 9 - 363x 12 + 141x 15 + x 18 ), (A.7) 
•fa = ^ 4 (-l - 104x 3 - 952x 6 - 4124x 9 - 3430x 12 

-1544x 15 - 7336x 18 + 616x 21 + 2x 24 ), (A.8) 

^ g = 3^ 3 (l _3 x 2 +x 3^ 1+3x 2 + 2x 3 + 9x 4 + 3a ,5 +:r 6^ 

x (1 + 9x 2 + 3x 3 + 18x 5 - 24x 6 + 9x 8 + x 9 ) 
x (1 - 9x 2 + 6x 3 + 81x 4 - 45x 5 - 39x 6 + 324x 7 
+153x 8 - 142x 9 + 486x 10 + 396X 11 + 582x 12 

+324x 13 + 198x 14 - 48x 15 + 81x 16 - 9x 17 + x 18 ), (A.9) 

^10 = -^2^5(1 - 177x 3 - 474x 6 - 7070x 9 - 104805x 12 

-542232x 15 - 862941x 18 - 1404072x 21 - 368055x 24 

+29 3 80x 27 - 55 2 84x 30 + 1173x 33 + x 36 ), (A.10) 
^ n = -l- 242x 3 + 605x 6 + 102729x 9 + 2270301x 12 
+17393277x 15 + 59389374x 18 + 189881835x 21 
+1106263389x 24 + 4869514969x 27 + 10595519759x 30 
+8054721004x 33 - 22319781x 36 - 4760052033x 39 
-8579472693x 42 - 1596123771x 45 + 66133914x 48 

-62045313x 51 - 1153603x 54 + 23221x 57 + llx 60 , (A.ll) 
^12 = V3^4(-2 + x 3 )(l-3x + 3x 2 + x 3 ) 

x (1 + 3x + 6x 2 + llx 3 + 12x 4 - 3x 5 + x 6 ) 

x (-2 - 32x 3 - 84x 6 - 134x 9 + x 12 ) 

x (1 + 6x + 12x 2 + 4x 3 + 45x 4 + 36x 5 + 60x 6 

+72x 7 - 45x 8 + 58x 9 - 48x 10 + 12x n + x 12 ) 
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x (1 - 6x + 24x 2 - 64a; 3 + 75x 4 + 456x 5 - 620x 6 
+252x 7 + 2070x 8 - 1618x 9 - 3072a; 10 + 3216a; 11 
+4003x 12 - 9696a; 13 + 1416x 14 + 11396x 15 
+1548x 16 - 5058x 17 + 460x 18 + 1632x 19 

+1653x 20 + 692x 21 + 192x 22 - 12x 23 + x 24 ), (A.12) 
= (i + i6x 3 + 96x 6 + 13x 9 + 13x 12 ) 

x (1 - 354x 3 - 17247x 6 + 92420x 9 - 6264417x 12 
-91630974x 15 - 414038735x 18 - 631690011x 21 
+3596512338x 24 + 43118516972x 27 + 215967505719x 30 
+533661527514x 33 + 582732421153x 36 + 284118813696x 39 
+450924775284x 42 + 1313707269872x 45 + 1846766455056x 48 
+403474854555x 51 - 263110973327x 54 - 22534762701x 57 
+685417938x 60 - 111537892x 63 - 798438x 66 

+57 48x 69 + x 72 ), (A. 13) 

^» 14 = foifoQ - 48x 3 - 741x 6 - 1924x 9 - 363x 12 + 141x 15 + x 18 ) 
x (1 + 504x 3 + 2421x 6 + 5676x 9 + 166356x 12 
+3098475x 15 + 22597638x 18 + 56826270x 21 
-73281168x 24 - 582904249x 27 - 862862121x 30 
+133470252x 33 + 317907519x 36 - 632536713x 39 
-77646699x 42 - 41502855x 45 - 2997252x 48 

+88 47x 51 + x 54 ), (A. 14) 

= V3^5(-5 + 65x 3 + 685x 6 + 3410x 9 + 11425x 12 

+57 35x 15 + 3 145x 18 - 520x 21 + x 24 ) 
x (1 - 6x + 6x 2 + 44x 3 + 21x 4 - 21x 5 + 676x 6 

+9x 7 - 9x 8 + 569x 9 + 2841x 10 - 2841X 11 

-1694x 12 + 13119x 13 - 13119x 14 + 10019x 15 

-4284x 16 + 4284x 17 + 4591x 18 - 1446x 19 + 1446x 20 

-496x 21 - 24x 22 + 24x 23 + x 24 ) 
x (1 + 6x + 30x 2 + 124x 3 + 279x 4 - 495x 5 + 3036x 6 

+2871x 7 - 2790x 8 + 60959x 9 - 13686x 10 - 19695X 11 

+469946x 12 - 200034x 13 + 128295x 14 + 602229x 15 

-2440926x 16 + 3056445x 17 - 422129x 18 - 9809094x 19 

+18607485x 20 + 20165779x 21 + 8262864x 22 

+74286585x 23 + 94839246x 24 + 71549460x 25 

+150594579x 26 + 118349119x 27 - 3156510x 28 

+30275751x 29 - 36357239x 30 - 138954870x 31 

-1389186x 32 + 73952184x 33 + 20894985x 34 

-28859229x 35 + 22894661x 36 + 16500675x 37 

-244451 lx 38 - 2237686x 39 + 1693800x 40 

+672156x 41 + 324606x 42 + 58950x 43 + 11034x 44 

-416x 45 + 600x 46 - 24x 47 + x 48 ). (A. 15) 
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A. 2 y 2 = x(x 2 - 1) 



V>i = 1, 

^2 = -2y, 

^ 3 = 3(-2 + x)x 2 (2 + x), 

^/, 4 = -Ay x 2 (-6 - 15x 2 + x 4 ), 

^ 5 = x 4 (-192 + 1632x 2 - 496x 4 



220x 6 



+ 5x 8 ), 



(A. 
(A 
(A 
(A 
(A 



16) 
17) 
18) 
19) 
20) 



A = -6y (-2 + x)x 6 (2 + x) 

x (-336 + 912x 2 - 1348x 4 - 100x 6 + x s ), (A.21) 

^7 = x 8 (27648 + 483840x 2 - 2951424x 4 + 2595456x 6 - 1101888x 8 

+447840x 10 - 31376x 12 - 1544x 14 + 7x 16 ), (A.22) 

^ 8 = ~8y x 10 (-6- 15x 2 + x 4 ) 

x (-18432 + 603648x 2 - 2432640x 4 + 2577312x 6 - 

702392x 8 + 47744x 10 - 23070x 12 - 412x 14 + x 16 ), (A.23) 

ifo = -3(-2 + x)x 14 (2 + x) 

x (16367616 - 154607616x 2 + 1527054336x 4 - 5301780480x 6 
+4162000896x 8 + 567207936x 10 - 1938695936x 12 + 731321472x 14 

-1489472x 16 + 5367072x 18 - 164000x 20 - 2316x 22 + 3x 24 ). (A.24) 

A.3 y 2 = x 2 (x + 1/4) 

V>i = 1, (A.25) 

4>2 = -2y, (A.26) 

^ 3 = x 3 (l + 3x), (A.27) 

ip 4 = -2yx 5 (l + 2x), (A.28) 

= x 10 (l + 5x + 5x 2 ), (A.29) 

ip 6 = -2y x 14 (l + x)(l + 3x), (A.30) 

fa = x 21 (l + 7x + 14x 2 + 7x 3 ), (A.31) 

%!>% = -2yx 27 (l + 2x)(l + 4x + 2x 2 ), (A.32) 

fa = x 36 (l + 3x)(l + 6x + 9x 2 + 3x 3 ), (A.33) 

fa = -2y x 44 (l + 3x + x 2 )(l + 5x + 5x 2 ), (A.34) 

fa x = x 55 (l + llx + 44x 2 + 77x 3 + 55x 4 + llx 5 ), (A.35) 

fa 2 = -2y x 65 (l + x)(l + 2x)(l + 3x)(l + 4x + x 2 ), (A.36) 

fa 3 = x 78 (l + 13x + 65x 2 + 156x 3 + 182x 4 + 91x 5 + 13x 6 ), (A.37) 

fa 4 = -2y x 90 (l + 5x + 6x 2 + x 3 )(l + 7x + 14x 2 + 7x 3 ), (A.38) 

fa 5 = x 105 (i + 3 x )(i + 5 x + 5 x 2 )(i + 7 x + i4x 2 + 8x 3 + x 4 ) (A.39) 

ip w = -2y x 119 (l + 2x)(l + 4x + 2x 2 )(l + 8x + 20x 2 + 16x 3 + 2x 4 ). (A.40) 
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